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Within the framework of effective action QED, we derive 
the light cone condition for homogeneous non-trivial QED 
vacua in the geometric optics approximation. Our result gen- 
eralizes the "unified formula" suggested by Latorre, Pascual 
and Tarrach and allows for the calculation of velocity shifts 
and refractive indices for soft photons travelling through these 
vacua. Furthermore, we clarify the connection between the 
light velocity shift and the scale anomaly. This study mo- 
tivates the introduction of a so-called effective action charge 
that characterizes the velocity modifying properties of the 
vacuum. Several applications are given concerning vacuum 
modifications caused by, e.g., strong fields, Casimir systems 
and high temperature. 

12.20.-m, 41.20.Jb, ll.10.Wx 



I. INTRODUCTION 

The vacuum considered as a medium has become 
a popular picture in quantum field theory. With re- 
servations due to the lack of understanding of non- 
perturbative vacuum phenomena, it is astonishing that 
analogies between the quantum vacuum and classical me- 
dia are frequently useful. 

One particular example is represented by the propa- 
gation of light in a vacuum which is modified by var- 
ious external environments, e.g., electromagnetic (EM) 
fields, temperature, geometric boundary configurations, 
gravitational background and non-trivial topologies. The 
concept of drawing the analogy is common to all of these 
cases: vacuum polarization allows the photon to exist as a 
virtual e + e _ -pair on which the various vacuum modifica- 
tions can act. Under certain assumptions, this influence 
on the loop process can effectively be described by an 
immediate influence of a (generally non-linear) medium 
on the photon itself, e.g., by refractive indices. This pro- 
gram was carried out among others by Adler Q , Brezin 
and Itzykson Q for magnetic fields, by Drummond and 
Hathrell 0] for gravitation, and by Scharnhorst Q] and 
Barton |5j] for a Casimir configuration. Further impor- 
tant examples are found in refs. 



A new physical insight into the phenomenon of photon 
propagation in non-trivial vacua has been given by La- 
torre, Pascual an Tarrach Comparing the known ve- 
locity shifts arising from different vacuum modifications, 
they were able to identify an intriguing general, so-called 
"unified" formula covering all these casesQ. They con- 
cluded that the polarization and direction averaged ve- 
locity shift is related to the (renormalized) background 
energy density u with a "universal" numerical coefficient 
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a 
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where m denotes the electron mass and a ~ 1/137. (In 
the case of gravitation, one a has to be replaced by 
the combination (G^m 2 ) involving Newton's constant). 
However, a complete derivation of the "unified formula" 
has not been given up to now. 

In the case of gravitation, light was shed on the prob- 
lem by Shore [lC| who proved a polarization sum rule that 
represents a generalization of eq.(|l]). Furthermore, he 
pointed out that the "universal" coefficient in eq.(|l|) can 
be related to the trace anomaly of the energy- momentum 
tensor in the case of weak EM background fields. 

One of the most remarkable features concerning vac- 
uum induced velocity shifts certainly is the fact that 
Sv > is not intrinsically forbidden in quantum field 
theories. This seems to offer the possibility of superlu- 
minal propagation, e.g., in curved spaces and Casimir 
vacua. Both examples share the property of a possible 
negative energy density u in eq.(^). 

The two questions, whether the signal (=wave front) 
velocity indeed exceeds c and whether superluminal prop- 
agation is observable in principle, could be resolved by 
calculating the velocity shift in the infinite frequency 
limit. But this is presently out of reach, because a resum- 
mation of the derivative expansion has to be achieved. 
However, without being able to answer these questions, 
let us just say that we find no grounds for violation of 
(micro-) causality in accordance with [p|P, p4]| , pl| . For a 
causality violation, a space-like signal and Lorentz in- 
variance (in the gravitational case: strong principle of 
equivalence) are necessary conditions. The latter is ex- 
plicitly violated in the above-mentioned examples. For 



1 In fact, the results of |s| cannot be embedded in the "uni- 
* E-mail address: holger.gies@uni-tuebingen.de ned formula". The solution to this problem is an aim of the 

present work. 
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an excellent discussion, the reader is referred to the work 
of Shore @. 

In the present work, we confine ourselves to the case of 
non-trivial vacua modified by QED phenomena. Within 
the effective action approach Jl2| , we derive a covariant 
light cone condition in sect. II which turns out to be a 
generalization of the "unified formula" . The necessary 
assumptions are analysed in detail. 

In this framework, we are able to clarify the relation 
between the velocity shift and the trace anomaly in sect. 

III. Our findings do not unveil a natural and physically 
meaningful connection. An alternative physical picture 
of the "universal" pre-factor is given instead which is 
called: effective action charge. Several applications of 
our light cone condition concerning EM fields, Casimir 
configurations and temperature are elaborated on in sect. 

IV. In the low-energy domain, we can easily recover all 
of the well-known results described by the "unified for- 
mula". However, the "universal constant" turns out to 
be neither constant nor universal when we drop the low- 
energy restriction. Instead, the concept of an effective 
action charge provides for an intuitive understanding of 
the velocity shifts at arbitrary energies. 

Conclusions are drawn in sect. V. 



II. THE LIGHT CONE CONDITION 

Consider light propagation in a non-trivial QED vac- 
uum (we will specify this terminology soon) characterized 
by a certain energy scale. Suppose that there exists an 
effective action which takes into account any QED quan- 
tum phenomena on higher scales and hence provides for 
an exact description of the propagation. In principle, this 
effective action will depend on any gauge and Lorentz in- 
variant scalar which we can construct. Throughout the 
paper, we will stick to the following essential assump- 
tions: 

1) The propagating photons characterized by are 
considered to be soft. This is equivalent to calculating 
the properties of the vacuum in the limit uj/m <C 1 where 
the scale is set by the Compton wavelength. 

2) The vacuum modification is homogeneous in space and 
time (but not necessarily isotropic). 

Referring to these assumptions, we can neglect any 
term in the effective action that involves derivatives of the 
field, since a derivative either acting on the background 
field vanishes (assumption 2)) or acting on the photon 
field f^ v contributes terms of the order 0{ui 2 /m 2 ) to the 
equation of motion. In the latter case, it is negligible 
because of assumption 1). 

We furthermore assume that 

3) Vacuum modifications caused by the propagating light 
itself are negligible. 

Assumption 3) justifies a linearization of the equations 
of motion with respect to f^ v but does not stand on the 
same footing as the former assumptions, since it is not 



essential for the formalism. Note that we do not demand 
that the deviation from the Maxwell Lagrangian should 
be small, corresponding to small vacuum modifications.^] 
Since it is unwieldy to establish a general formalism for 
arbitrary numbers of Lorentz vectors and tensors char- 
acterizing the vacuum, we first consider a vacuum only 
modified by EM fields. Hence, the dynamical building 
blocks of the effective action which respect Lorentz and 
gauge invariance are given by the field strength tensor 
and its dual 

F ^ = qh A v _ qv A v ( 2a ) 

The lowest-order linearly independent scalars are 

x:=\F^F^ = \{B 2 -E 2 ) (3a) 
y := *F^=E-B, (3b) 

The normalization is chosen in such a way that the 
Maxwell Lagrangian can be written Cm = — By tak- 
ing advantage of the antisymmetry of F 11 " and by virtue 
of the relations |Q 

Ft tap» a _ * F u^ = 2x g^, (4a) 
F^ a *F v a = W a = y g^ , (4b) 

using the metric g — diag(— , +, +, +), it is easy to verify 
(i) the vanishing of odd-order invariants and (ii) that in- 
variants of arbitrary order can be reduced to expressions 
only involving x n y m where n,m = 0,1,2.... Besides, 
note that parity invariance demands for m to be even. 

Consequently, the complete effective action becomes 
extremely simplified, turning out to be a function of x 
and y only. The corresponding Lagrangian reads 

C = C(x,y). (5) 

We obtain the equations of motion from C by variation 



11 d(d fl A^) dA» 
= d„(d x C F» v + d y C *F» V ) (6) 

where d x , d y denote the partial derivatives with respect 
to the field strength invariants (||) (and should not be 
confused with space-time derivatives 9 M ). 

If we take advantage of the Bianchi identity while mov- 
ing <9 M to the right, we arrive at 



2 In principle, the Lagrangian C can contain imaginary parts 
indicating the instability of the modified vacuum. In the fol- 
lowing, it is understood that we take into account only the real 
part of C which is solely responsible for the field equations. 

3 x and y are usually called T and Q. We do not follow this 
convention for reasons of simplicity. 
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(7) 



where M£ g is given by 



M M£ . = F M^ a/3 ( 9 2 £) + * Fq/3 (5 2 £) 

+&„£ (F^ *Fa3 + *F^F ctfi ) . (8) 

Note that M is antisymmetric in the upper as well as the 
lower indices: M% = -M% = M%. 

In general, contains background fields F^ u and the 
propagating photon field f^ u . According to assumption 
2), the derivative acting on F^ v vanishes 

d»F XK = 0^f XK . (9) 
Inserting eq.(^), eq.(^|) yields in Fourier space 

o = (d x c) k^r + (|m$) k^r? . (io) 

Introducing a gauge potential a M for the propagating field 
f^ v , we may write 

= k^a" - k v a^ = a (k"e u - jfe"e") , (11) 



where a := ^/a^a^ and e M = a^/a. Here, the polarization 
vectors e M are normalized to 1. 

Establishing the Lorentz gauge fc M e^ = 0, we get 



= (d x C) k 2 e» + M^k^k a eP, 



(12) 



where we used the antisymmetry of M£ g . 

The next important step is to multiply eq.(|T2|) by e v 
and average over polarization states according to the 
well-known rule 



9 
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(13) 



pol. 



where the additional terms on the right-hand side of Jl3] ) 
vanish with the aid of the antisymmetry of M^ g . We find 
for eq.@ 

= 2(8 X C) k 2 + MZ k^k a . (14) 

Equation ([l4]) already represents a light cone condition 
and actually indicates that the familiar k 2 = will in 
general not hold for arbitrary Lagrangians. Our final 
task is to put MZ in a convenient shape. Using the 
powerful relations (^), we obtain 

MZ = 2 [\F^F av {dl + d 2 )C + 5Z(yd xy £ - xd 2 C)} . 

(15) 

Introducing the Maxwell energy-momentum tensor 

T» a = F^F au -x6Z, (16) 

this leads to 



+5»(\x{dl-dl)C + yd xv C)} 



(17) 



However, the Maxwell energy-momentum tensor in gen- 
eral is devoid of any physical meaning, since we are sim- 
ply not dealing with the Maxwell Lagrangian. The right 
quantity to deal with is therefore the vacuum expectation 
value (VEV) of the energy-momentum tensor defined byn 



-gSg t 



- ,T:= / #Xy/=gC.. 

oguv J 



(18) 



where T denotes the effective action. Performing the cal- 
culation, we arrive at 

{T^)*y = -T^{d x C) + (C - xd x C - yd v C) . (19) 

Solving eq. ( p"9| ) for T^ u and inserting into eq. ( |l7j ) , we can 
present MZ in its final form 



d x C 



2 



(T^) xy + 5^x(d 2 x -d 2 )C 



Ud 2 x +d 2 )C 



(C-xd x C-yd y C) 



(20) 



Substituting MZ into eq.(^), we end up with the de- 
sired light cone condition for EM field modified vacua 
fulfilling the above-mentioned assumptions 



where 



Q = 



m+d 2 )c 



(21) 



(22) 

To extend the validity of the light cone condition to arbi- 
trary non-trivial vacua, we have to take the vacuum ex- 
pectation value of eq.([jl"l) with respect to the additional 
vacuum modifications parametrized by the (collective) 
label z 

k 2 = z(0\ Q (T^jxy |0) z k^ku . (23) 
Inserting a complete set of intermediate states, we obtain 

k 2 = Y, t (0\Q\i),MT' M ')xu\0).k^k v . (24) 



4 Note that the variation with respect to the metric tensor 
is just a trick to calculate the symmetric energy-momentum 
tensor. With some care, the same result can be obtained by 
canonical methods. 
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In the following, we consider the vacuum to behave as a 
passive medium in which the EM fields and the further 
vacuum modifications z remain in a state of static equi- 
librium. Since Q solely depends on x and y (via C(x, y)), 
this assumption of passivity leads to 



Q\i) z = (Q} z S Ql . 



(25) 



Equation ([2q) states that the vacuum exhibits no back- 
reaction caused by the EM fields while switching on z. 

Q depends functionally on L{x,y), which is, as usual, 
defined via the functional integral over the fluctuating 
fields. Taking the expectation value of Q hence leads 
back to integrating over the field configurations which 
respect the modified vacuum. E.g., if the modification z 
imposes boundary conditions on the fields, the functional 
integral has to be taken over the fields which fulfil these 
boundary conditions. Therefore, taking the VEV of Q 
defines the new effective Lagrangian characterizing the 
complete non-trivial vacuum 



(Q) z = (Q(C(x,y))) z =Q(C(x,y;z)) 



(26) 



We finally arrive at the light cone condition for arbitrary 
homogeneous non-trivial vacua 



k — Q{x, y, z) (T^ )xyz k^ku 



(27) 



Remember that the validity of the light cone condition 
eq. (|27]) is not restricted to results of perturbation theory 
or only small modifications of the Maxwell Lagrangian. 
It is an exact statement in the sense of effective theories. 

Now, the terminology "modified QED vacuum" should 
be clarified: from the derivation of the light cone con- 
dition, it is obvious that the implicit space-time depen- 
dence of C should only be contained in the field variables. 
Furthermore, the vacuum has to fulfil the demand for 
passivity. Otherwise the light cone condition (^) only 
represents a zeroth order approximation of the infinite 
sum over intermediate states in eq.(p5|). 

As a third remark, we want to point out that the sum 
over polarization states is not necessary for the derivation 
of a light cone condition. By summing, we even exclude 
the study of birefringence from the formalism which is 
certainly the most important experimental application 
jll] |l6|] . But for a projection on the polarization eigen- 
states, the y ra -terms have to be rewritten in terms of the 
field strength tensor which is practically impossible for 
arbitrary C. 

In the remainder of the section, we calculate further 
representations of eq. ( p7|) by choosing a certain reference 
frame and introducing 



Ifel 



k° . 
Ifel 



(v,k) , 



(28) 



where we defined the phase velocity by v := k°/\k\. For 
eq.(p7n, we obtain 



Equation ( p9| ) clearly demonstrates that the light cone 
condition is a generalization of the "unified formula" of 
Latorre, Pascual and Tarrach S. 

In general, the Q-factor will depend on all the variables 
and parameters of £ and hence will naturally be neither 
universal nor constant. Besides, the daunting structure of 
the Q-factor will simplify in the case of small corrections 
to Cm- As will be shown in sect. IV, the denominator 
then reduces to 1. 

Another representation of the light cone condition is 
found by averaging over propagation directions, i.e., in- 
tegrating over k € S 2 



i-Q(i(r°°) + i(r°)) 

1 + Q (TOO) 



(30) 



For Q(T 0Q ) <C 1 and (T a a ) being even of lower order, this 
reduces to 



v 2 = l-\Q{T m ) = l- A -Q Ul 



(31) 



v 2 = 1 - Q (T^k^h ■ 



(29) 



where u denotes the (renormalized) energy density of the 
modified vacuum. 



III. VELOCITY SHIFT AND SCALE ANOMALY 

In his paper, Shore suggested a deeper connection 
between the velocity shift and the scale anomaly. For 
the Heisenberg-Euler Lagrangian, he showed that the co- 
efficients of the x 2 and y 2 terms in the scale anomaly 
are precisely those appearing in the velocity shift for the 
different polarization states. 

Within the framework developed so far, we will at- 
tempt to clarify the relation between the scale anomaly 
and the velocity shift. Therefore, we have to investigate 
whether the terms in the Q-factor can be expressed in 
terms of the anomaly. For reasons of simplicity, we limit 
this consideration to the case of a purely EM field mod- 
ified vacuum. From eq . (|l9|) , we can read off the scale 
anomaly 

(T a a ) = 4(£ - xd x C - yd y C) . (32) 

By differentiation, we find 

dx (T a a ) = -4(xd 2 £ + yd xy £) , (33a) 
d y (T a a ) = -±(yd 2 £ + xd xy C) . (33b) 

From eqs.(p3|) immediately follows 

(d 2 x +d 2 )C = -(I + i) dxy C-\{ld x + l -d v ){T« a ) . 

(34) 

This expression is proportional to the numerator of the 
Q- factor eq. (p2[) . Using similar techniques, we can also 
rewrite the denominator, but the result is not very illu- 
minating 
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denom.(Q) = (d x C) 2 - (d x C){% {8 2 x + 8 2 y )C + \d x {T a a )) 



+\[(% + %)C]{T a a ) 



(35) 



Fortunately, approximating eq.(p5[) by 1 will be appro- 
priate to the applications of sect. IV. 

It is already obvious from cq.(|34|) that there is no 
immediate connection between (T a a ) and the velocity 
shift eq.(29). The findings of Shore arise from the spe- 
cial structure of the Heisenberg-Euler Lagrangian where 
d xy C = 0. In general, higher-order mixed terms are not 
forbidden by gauge, Lorentz or parity invariance. Re- 
ferring to eq.([M|), the introduction of the scale anomaly 
appears to be artificial rather than interpretable. Even 
Shore's conjecture that the sign of the scale anomaly is 
linked to the sign of the velocity shift cannot be main- 
tained. 

Instead, we favour the pure effective action formula- 
tion, i.e., the left-hand side of eq.(|34|), since it offers a 
new intuitive picture. Referring to eq.(^9|), the value and 
sign of the velocity shift result from the competition be- 
tween the VEV of the energy-momentum tensor and the 
Q-factor. Both are a priori neither positive nor bounded 
by symmetry principles. Let us restrict the following in- 
vestigation to the case of small corrections to the Maxwell 
Lagrangian, i.e., 



V 2 C = 2Q. 



(36) 



Due to the similarity to the (2-D) Poisson equation, we 
will call Q from now on the effective action charge in 
field space. The classical vacuum Cm = —x is uncharged 
and hence v = 1. As we will soon demonstrate, the 
pure QED vacuum has a small positive charge at the 
origin in field space (x = y = 0). For increasing field 
strength, {T<- iU ) certainly also increases without upper 
bound, so we expect Q to decrease in order to produce no 
unphysical velocity shift > 1. It is therefore reasonable 
to presume localized effective action charge distributions 
centred upon the origin in field space. The results of sect. 
IV will confirm this charge-like picture. 



IV. APPLICATIONS OF THE LIGHT CONE 
CONDITION 

Up to now, the light cone condition might be regarded 
as a nice frame without a picture enclosed, since it is 
much easier to talk about all-loop or non-perturbative 
effective actions than to calculate one. 

Indeed, the effective actions which we are going to in- 
sert will not reach beyond two-loop order. Their general 
structure can be characterized by 



C = £ M + C c ; < 1 , 

where C r contains the correction terms. 



(37) 



Regarding the denominator expression of the effective 
action charge ((3^), the scale anomaly (T a a ) is of the same 
order as C c . Hence, eq.(p9) simply reduces to 



eq.(35) = 1 + 0{C C ), 



(38) 



and the approximation Q = iV 2 £ = iV 2 £ c is justified. 



A. Weak EM Fields 

According to the authors of ref . jL7) , the two- loop cor- 
rected Heisenberg-Euler Lagrangian (weak-field limit of 
the complete one-loop approximated effective QED La- 
grangian) reads 



C = —x + ci x 2 + C2 y 2 ■ 



where 



Cl 



C-2 



8a 2 
45to 4 
I4a 2 
45m 4 



40 a 
1315 a 



252 7T 



(39) 

(40a) 
(40b) 



With the aid of the light cone condition eq.(pl|), we im- 
mediately obtain for the polarization and propagation 
direction averaged velocity (v = Vvr*) 



c\+c 2 , 
4a 2 



1 



135m 4 



11 



1955 a 



36 7T 



-(E 2 

2 V 



B 2 ) 



(41) 
(42) 



In the well-known one-loop part of eq. ( |42|) , we can iden- 



tify the factor of .tt a 4 as the "universal constant" of 



135/71 



the "unified formula" eq.(Q). At this stage, it is already 
understandable that all of the known QED induced ve- 
locity shifts share this universal factor in the low-energy 
limit, since they are all based on the Heisenberg-Euler 
Lagrangian. Even the results in gravitation involve the 
same e + e~-loop calculation (of course, in a curved space- 
time ||]).| It is furthermore obvious that the two-loop 
correction ^§|p^r is as universal as the number 11. (Note 
that modifications from the denominator of Q eq.(|35|) 
contribute to the order C(a 4 )). 



B. Strong Magnetic Fields 

Since the Heisenberg-Euler Lagrangian eq.(^) repre- 
sents a weak field limit, eq.([i"l|) denotes the value of 
the effective action charge at the origin in field space 



5 In particular, there is nothing mysterious about the factor 
11 as it is sometimes found in the literature. 
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(x, y = 0). In this subsection, we analyse the form of Q 
along the positive ir-axis (pure magnetic fields). As our 
starting point, we use Schwinger's famous formula for the 
one-loop effective QED Lagrangian |l3[ 



lOO 

1 fds 



8tt 2 j s 
o 



(es) 2 \y\ coth\es(\/ x 2 +y 2 + x) 

/ is 2 

x cot[es(y x 2 +y 2 — x)?J — -(es) 2 x 



(43) 



It is understood that the convergence is implicitly en- 



sured by the prescription m 



ie. (Note that we 



have not performed a proper time Wick rotation yet.) 

It will be useful to reparametrize the field space with 
new coordinates 



a := (\/ x 2 + y 2 + x) 2 , b := (\/x 2 + y 2 — x) 2 
\y\=ab 1 x= 1 1 (a 2 -b 2 ). 

The Laplacian in terms of a and b reads 

_o 1 



-b 2 



{dl + d 2 b 



(44a) 
(44b) 



(45) 



For the term in the square brackets in eq.(|43]), we easily 
find 



a 2 + b 2 
2(es) 2 
a 2 + b 2 



(d 2 + d 2 ) [a&coth esa cot esb] 



sb cot esb 



sinh esa 



(esa coth esa — 1) (46) 



esa coth esa 

H ; 5 ( es " C °t — 1 ) 

sinh esb 



Q(h) 



1 a 



(2/i 2 - + - h - 3/i 2 - 4ftliir(/i) 

(49) 



-2/iln27r+i + 4C'(-l,/i) + -^ 



where V 1 denotes the logarithmic derivative of the T- 
function and £' is the first derivative of the Hurwitz Zeta 
function with respect to the first argument p3fl . 

For strong fields, the last term of eq.(^9|) oc ^ oc \B\ 
dominates the expression in the square brackets. Hence, 
the effective action charge decreases with 



Q(B) 



la 1 1 



, for B — > oo 



6 7r B CT B 

which supports the charge picture (fig. 1). 



(50) 
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FIG. 1. Effective action charge Q{B) — in units of 
versus magnetic field B in units of the critical field strength 

m 2 



Confining ourselves to purely magnetic fields (x — \B 2 , 
y = 0=>6 = 0, a — \B\), we obtain 



eq.(46) 



2(es) 



2 r 



esa coth esa — 1 I 
sinh 2 esa 3 



esa coth esa 



(47) 



The complete formula for the effective action charge 

2 

might be written (substitution: z := esa, h := = 



2B > 



loo 

1 a f dz _ 2hz 



z coth z — 1 1 

— z coth z 
3 



sinh 2 z 



(48) 



With some effort, the evaluation of the integral can be 
performed analytically by standard means of dimensional 
regularization. Details are given in App. A. The result 
is 



The contraction of the energy-momentum tensor VEV 
may be cast into the form 

(T^kftkv = B 2 (B- kf + O(a) = B 2 sin 2 + O{o) , 

(51) 

where 9 measures the angle between the B-field and the 
propagation direction. 

Finally, the light cone condition cq.(|29|) yields for ar- 
bitrary background fields consistent with the one-loop 
approximation (h = -fgr) 



= 1 - 



a sin 2 6 



7T 2 

x 



\2B 2 3jV [ - 1+ 2b) B ml \2B) 4S 2 



B cr 
2B 



%ln27r+-+4C'(-l,f») + 



g ■ -s \ -i 2B > ~ 3B cr 



(52) 



G 



The first derivative of the Hurwitz Zeta function at —1 
can be related to the generalized T-function of first order 

ri M 



C(-l,h) = -hlnh + lnTi(l + h) - L x + 



1 

12 



(53) 



where L\ = 0.248 754477. . . is a pure number and can 
be obtained from the Raabe integral |lij . 

Using eq.(|53|), one can show that eq.(|52]) is identical 
to the findings of Tsai and Erber Equation (52) is 
plotted in fig. 2. Although the velocity shift increases 
proportional to the magnetic field for large B, the total 
amount of the velocity shift remains comparably small 
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FIG. 2. Square velocity v 2 versus magnetic field B in units 

2 

of the critical field strength B CT — — . The dashed curve 
indicates the region where two-loop corrections become im- 
portant. 



5v ~ 9.58.. • 10" 



at 



B = B CI = — (54) 
e 



for strong B-fields consistent with the one-loop approx- 
imation, i.e., -jp- < — ~ 430. Taking higher-order 
loop calculations into account, we expect a stronger de- 
crease of Q(B) for large B in order to let QiT^k^k^ be 
bounded. 



C. Casimir Vacua (Scharnhorst effect) 

One curious result regarding vacuum induced velocity 
shifts is the possibility of superluminal phase and group 
velocities. As mentioned above, e.g., Casimir vacua can 
create positive velocity shifts, since a negative shift of the 
zero point energy is permitted. For the configuration of 
perfectly conducting parallel plates of distance a, (T^") 
is found to be 20 



(T 



liv\ 



720a 4 



(55) 



where the symmetry axis points along the 3-direction. 
The effective action charge has to be evaluated in the 
zero-field limit. In concordance with experimental facil- 
ities, the plate separation a is treated as a macroscopic 
parameter (a oc fJ-m); otherwise, we would violate the 
soft-photon approximation, since the photon wavelength 
A has to obey A <C a to validate the concept of treating 
the Casimir region as a (macroscopic) medium. 

The magnitude of a implies that we can neglect the 
a-dependence of Q which is exponentially suppressed by 
ma 1 (this point will become clearer in the following 
section). Here, Q is simply given by eq.(41) 



ci + c 2 



2a 2 
45m 4 



11 



1955 a 



36 7T 



from which directly follows using eqs.(55) and (£ 
la 2 / 1955 a\ it 2 



v = l 



(90) 2 to 4 



11 + 



36 7r / a 4 



(56) 



(57) 



for a propagation perpendicular to the plates (a parallel 
propagation will, of course, not be modified). 

Equation ([57]) represents the two-loop corrected ver- 
sion of Scharnhorst 's formula Note that the two- 
loop correction enhances the velocity shift. As was re- 
cently found by Kong and Ravndal pj| , the radiative 
correction to the Casimir energy is of order a 4 s : 



(T 



(Hl\ 



11 



7r 4 a 2 



720a 4 (90) 2 -30-16 



(58) 



At the two-loop level of eq. (p7|) , this correction can ob- 
viously be neglected. Even three-loop contributions in Q 
would be more important. But it is interesting to note 
that this correction also contributes positively to v. 



D. Finite Temperature 

In the remaining sections, we reveal the manifold fea- 
tures of temperature induced velocity shifts. Unlike the 
Scharnhorst effect, we do not recognize a principal ob- 
stacle against measurability here, and the results allow 
for an immediate physical interpretation. The following 
calculations are restricted to the one-loop level. 

We begin with the one-loop correction to the effective 
QED Lagrangian at finite temperature which can be de- 
composed according to 

C c (x, y, T) = C c (x, y, T = 0) + A£(x, y, T) , (59) 

whereby C c {x, y, T = 0) denotes the usual zero-tempera- 
ture Lagrangian eq. (|43|) . 

For purely magnetic fields, A£(x,y,T) was calculated 
by Dittrich p3j 



7 



lOO 

A£(B,T) = f ^e' m2s esB cot esB (60) 



S2 



xT 



e 2 (0,47risT 2 ) - 



1 



2Tw¥s 



The Jacobi 0-function is defined by |p3| 

OO 

9 2 (0,-g)= ^ exp(-ig(n+i)*) . ((ill 



n— — oo 



The effective action charge can be decomposed similarly 
to eq.(|59|) into 



(62) 



Q(x, y, T) = Q(x, y, T = 0) + AQ(x, y, T) 



= \ V 2 £ c (.t, y, T = 0) + | V 2 AC(x, y, T) . 

Q(x,y,T = 0) clearly corresponds to the zero-tempera- 
ture case as treated above. 

Since we have to differentiate with respect to x and 
y, it is not sufficient for the calculation of AQ to con- 
sider magnetic fields only in eq. (|60|) . Not until we have 
carried out the Laplacian are we allowed to set E = 0. 
Indeed, we have to take this limit E — ► in the end, be- 
cause the principle of equilibrium thermodynamics would 
otherwise be violated. Besides, the above-mentioned as- 
sumption of passivity of the vacuum is only fulfilled for 
E = 0. The appropriate expression is simply obtained 
by replacing 

esB cot esB (63a) 
by the gauge and Lorentz invariant terms 



(esa)(esb) coth(esa) cot(esfr) 



(63b) 



in analogy to eq.(f43|). Again, we make use of the advan- 
tageous coordinates a, b in field space defined in eqs.(44). 
The result of the differentiation was already found in 
eq.(^); hence we obtain for the temperature induced 
effective action charge for purely magnetic fields (a = 
5,6 = 0) 



AQ(B,T) 



lOO 

/7F a f ds 



esacothesa— 1 esa 
sinh 2 esa 3 

cT(e 2 (0,4.i S T 2 )-^) 



coth esa 



a 1 [ds 



it a j s 




esacothesa — 1 esa 



coth esa 



sinh esa 

00 2 
x^V i,m e-^. (64) 

n=l 

In the last line, we made use of the identity p2| 



9 2 (0,47risT 2 ) 
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(65) 



Our task is to evaluate eq.(|64|) in the various limits. First, 
we consider pure temperature phenomena with vanishing 
field strength. The temperature dependent part of the 
effective action charge reduces to 



22 



AQ(B = 0,T)= 2^a 2 / dsse~ 



22 a 2 



4 r )/n 4£(-l) n (^) 2 ^(». (66) 

71=1 

whereby we have taken advantage of the representation 

00 

2(f)" K v (ii) = J duu^exp (-«-£) (67) 


for the modified Bessel function and have rotated the 
contour. 

For low temperature, we may use the asymptotic ex- 
pansion of K2 (x) for i> 1 



K 2 {x) -- 
In this limit, we find 
AQ{B = 0,T^ 0) 



-(1 + 0®). 



(68) 



22 a 2 
45^ 



7r /m\ 2 



e T 



(69) 



Hence, the effective action charge is perfectly described 
by Q(B = 0, T = 0) = ci + c 2 eq.© in this limit, while 
the influence of temperature on Q vanishes as it should. 
(Note that in the case of Scharnhorst's effect a similar 
term AQ(B = 0,ma ^> 1) also vanishes by drawing the 
analogy T cx ~.) 

Next, we investigate the high-temperature limit 
T/m S> 1 of eq. (|66|) . The calculation is, however, much 
more involved, so we simply state the result 



22 a 2 

AQ(T > m) = —-7Z—7 

45 to 4 



4 r 4 V t 6 



(70) 

=const.. The interested 



where ki = 0.123 749 077 470. 
reader is referred to App. B. 

Therefore, we arrive at the remarkable result that the 
complete effective action charge 



Q(T > to) = Q{T = 0) + AQ(T > to) 

„2 



11 , Of ,„ / 



(71) 
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decreases rapidly oc 1/T 4 . Q(B = 0,T) is plotted in 
fig. 3. The influence of temperature causes the effective 
action charge to evaporate. Numerical results astonish- 
ingly indicate that eq. (|7l| ) is already a reasonable approx- 
imation for ~ 1.4 (error< 5%) where real e + e~-pair 
creation is energetically impossible and the vacuum is 
essentially modified by a photon gas. This excludes the 
interpretation that eq.(|7l|) is a pure threshold effect of 
pair production. 




10.0 



m/T 



FIG. 3. Effective action charge QiT) = Q(B = 0,T = 0) 
+AQ(T) in units of ff-^; for ni S h temperature, Q(T) de- 
creases proportional to =y. 



To complete the high temperature/zero field analysis of 
the light cone condition we need the VEV of the energy- 
momentum tensor which is given by (see, e.g., p4|) 

(T^)t = ^ (N B + lN F ) T 4 diag(3, 1, 1, 1) . (72) 

The integer variables Nb and N-p denote the number of 
bosonic and fermionic degrees of freedom at a given tem- 
perature. For QED, we obtain 

N B = 2 , N F = (73a) 

for T <C to (photon gas) 
-/Vb = 2 , N F = 4 (73b) 

for T>ra (photon + ultrarelativistic 
e + and e~ fermion gas) 

It is appropriate to employ eq. ([H]) for the light con e co n- 
dition. Using our findings in eqs.[(|69|), fl4l|), (f72]), (73a)], 
we recover the well-known result [|3| for low temperature 



44-7T 



2 t^4 
2 J 



2025 



(74) 



which according to fig. 3 is valid for Tim < 0.16 (error < 
5%) (T < 10 9 K). Substituting eq.flnj) into eq.Q and 
using eqs.[(|72])) d73b[ )] for T » to, we finally arrive at 



the velocity of soft photons moving in a photon and ul- 
trarelativistic e + e~ gas 



121 22 

V = 1 Kl 7T a 

8100 

= 1 - 9.72.. • 10~ 7 - 



0{^\ = const. 



(f^) ■ (75) 



In eqs.[([74|), (|75|)] , we found that the velocity shift in- 
creases proportional to T 4 for low temperature but ap- 
proaches a constant value in the high-temperature limit. 
This can be understood in terms of the effective ac- 
tion charge which evaporates sufficiently fast compared 
to the increase of the energy-momentum tensor VEV. 
Obviously, the shift described by eq.([75|) remains small; 
therefore the deviation from the vacuum velocity does not 
become seriously important (e.g., for the construction of 
cosmological models). However, one should keep in mind 
that, if the temperature exceeds the masses of further 
charged particles, each particle will contribute additively 
to Q and will increase the respective number of degrees 
of freedom Vb or Vp. 



E. Casimir Vacua at Finite Temperature 

The combination of thermal and Casimir phenomena 
is in itself worthwhile studying, because both effects en- 
ter the formalism via boundary conditions but lead to 
opposite results. In the following, we want to investi- 
gate where and why the respective effect dominates the 
velocity shift. The determining order parameter is the 
dimensionless combination Ta. Nevertheless, the plate 
separation a has to be considered as a macroscopic quan- 
tity (a ~ \irn) . 

First, we consider the low-temperature region. Accord- 
ing to Brown and Maclay |20|, the VEV of depending 
on a and T is given by 



irp00\a _ ^ 1 
V ,T 720 a 4 



(T 33 )£ 



240 a 4 



C(3) T 3 



for Ta -> . 



(76) 



for the parallel plate configuration (C(3) = 1.202056 . . .). 
The light cone condition eq. (|29|) for a propagation per- 
pendicular to the plates (fc M = (v, 0, 0, 1)) yields 



v = 1 



(90) 2 to 4 
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180C(3) 



(Taf 



(77) 



In the low-T limit, the (Ta) 3 -term can be neglected and 
we only rediscover Scharnhorst's result. But we do not 
find an additional velocity shift proportional to T 4 which 
could have been expected from eq.(|74|). This clearly 
arises from the fact that none of the (quantized) per- 
pendicular modes can be excited at low temperature. 







The (Ta) 3 -term in eg. (J77|) will become important for 
Ta = 0(1), i.e., T > 2000K for a ~ fim. This shows 
that the Scharnhorst effect is not sensitive to tempera- 
ture perturbations. 

For increasing temperature, we encounter an interme- 
diate temperature region characterized by the condition 
1 <C Ta -C ma which corresponds to 0.2eV^cC T <C 
0.5MeV. This implies that Q = Q(T = 0) is a justified 
approximation and the thermal contribution of an e + e~ 
gas does not have to be taken into account. 

Using further results of Brown and Maclay p0| 



(T° 



_rp4 
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(78) 
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for Ta > 1 , 



we find 



v = l 



Air 2 a 2 
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(45) 2 m 4 v " ' 36 w " ^ ( Ta ) 3 



45C(3) 



(79) 

In this limit, only the modifications caused by the black 
body radiation become important. A term proportional 
to 1/a 4 as a consequence of certain missing zero-point 
fluctuations does not occur, since higher (perpendicular) 
modes have been thermally excited. 

For T > m, we will certainly recover eq. (|75|) with neg- 
ligible l/(Ta) 3 Casimir corrections. Anyway, the concept 
of solid plates is (at least experimentally) meaningless in 
this domain. 



square brackets for z 3> 
manipulations of eqs.[(|66|) 



-z. Following the 



we arrive at 



AQ(T, B > B CI ) 



la 1 B 
ZnlPBT 



E(-i)"^ 



n=l 



(81) 



Note that it was not necessary to impose any conditions 
on the magnitude of T to arrive at eq. j8l|) . But, as men- 
tioned above, the field-temperature phenomena decouple 
for T <C m due to the asymptotic behaviour of K\(^n) oc 
exp(-f ) in eq.@; hence, AQ(f > 1,B > B CT ) -> 0. 
Using similar techniques as applied in App. B, the high- 
temperature limit of eq.(|8l|) can be determined. The 
result for T 3> m and B ^> B cl - is 



AQ(T, B) 
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(82) 



where fc 2 = 0.213 139 199 408 . . . =const. (see eq.(|B14|). 
To obtain the complete effective action charge Q, we add 
the strong-field contribution Q(T — 0) which was found 
in eq.(|50|) 



Q(T > to, B » B c 



la e 
6n 2 BT^^ 

2 a 2 1 
— ko — 7 - - 
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(83) 



F. Finite Temperature and Magnetic Fields 

For low temperature as well as for weak fields, thermal 
phenomena decouple from magnetic vacuum modifica- 
tions, because the effective action charge is not sensitive 
to weak influences. The velocity shifts can simply be de- 
scribed by an addition of the respective above-calculated 
ones. The only non-trivial interplay can be found in the 
domain of strong fields in hot surroundings (e.g., hot neu- 
tron stars). Our intension is to evaluate the thermal ef- 
fective action charge contribution given in eq.dfMf) in this 



limit. Therefore, we substitute z = esa (h — ^— = -fg-) 

3 v lea AB ' 



where we have introduced the convenient dimcnsionlcss 
variables B = -M- = ^1 and T — — which satisfy 

73 cr m J 

B,T ^> 1. Equation (l83|) describes the same features of 
the effective action charge which we have encountered in 
previous examples: Q is centred upon the origin in field 
space, decreases proportional to 1/B and evaporates with 
increasing temperature. 

To calculate the velocity shift, we need the energy den- 
sity which consists of three parts 
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cq.(75) 



AQ(h, T) 
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(80) 



In this representation, it is obvious that the integral is 
dominated by small values of z for weak fields (h 1) 
and vice versa, i.e., large z for strong fields. We are 
interested in the latter, so we expand the term in the 



The last term of eq.(|8J) is connected with the Lagrangian 
via the free energy (density) according to 



A(T UU )£ = F + TS = F -T 
r)C B 

- C B I T T 



dF 
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The leading mixed contribution C B to L is found in ref. 
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A(T 00 )f = C JLt 2 . (86) 



We finally arrive at the polarization and propagation di- 
rection averaged velocity shift for strong fields and high 
temperature 
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(87) 
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At T 2 / B = 1.74 . . ., we find a minimal velocity shift 

\Sv\ ~ 3.20- 1CT 5 . (88) 

At the same time, this number approximately sets the 
scale of a typical velocity shift for strong fields consistent 
with the one-loop a jproximation. This is also confirmed 



by the result of eq. 



V. CONCLUSIONS 

In this work, we studied light propagation in non- 
trivial QED vacua in the geometric optics approxima- 
tion. For any given QED effective action describing soft 
photons, we derived the light cone condition averaged 
over polarization states. This result generalizes the "uni- 
fied formula" found by Latorre, Pascual and Tarrach j|] 
which turned out to be the low-energy limit of our light 
cone condition. 

We furthermore clarified the connection between light 
velocity shifts and the scale anomaly suggested by Shore 
|jTof . Unfortunately, our findings do not indicate an im- 
mediate connection hinting at deeper physical grounds. 

Instead, the structure of the light cone condition sug- 
gests introducing the intuitive physical picture of an ef- 
fective action charge Q showing a localized profile in field 
space centred upon the origin. This charge directly char- 
acterizes the properties of the modified vacuum which 
are responsible for velocity shifts. 

Within this conceptual framework, we analysed sev- 
eral modified QED vacua and calculated the respective 
modified velocities. The inverse velocities are equal to 
the refractive indices of the modified vacua in the low- 
frequency limit. Hence, these velocities are phase as well 
as group velocities - the latter due to their independence 
of frequency. In the low-energy limit, we recovered all 
known results which were already perfectly described by 
the "unified formula". 

For arbitrary magnetic fields, we reproduced the find- 
ings of Tsai and Erber Q using our comparably simple 
formalism. 

In the sequel, we calculated the next-to-leading order 
corrections to the Scharnhorst effect. 

Finally, we investigated the influence of temperature 
on the velocity shifts. The evaporation of the effective 



action charge turns out to be the dominating effect in 
the high-temperature domain. It causes the velocity shift 
to approach a constant value. Only when a strong mag- 
netic field is involved, does the light cone condition fail 
to provide for a bound of the velocity shift in our exam- 
ples. But we expect higher-order loop corrections to stop 
an unbounded growth of the velocity shift by inducing a 
faster decrease of the effective action charge far from the 
origin in field space. 

Referring to the light cone condition, the sign of the 
velocity shift is in general determined by the sign of the 
effective action charge and the vacuum energy density. 
However, up to now, we have not been able to construct 
an example which exhibits a negative effective action 
charge in QED. This might be a general characteristic 
of the abelian theory. Indeed, the one-loop effective ac- 
tion of a covariant constant chromomagnetic background 
field 1 25 1 (naively) possesses a negative effective action 
charge. 

We would like to conclude with the remarkable obser- 
vation that parity violating terms in the effective action 
proportional to y 2n+1 , n = 0, 1, 2 . . . will not contribute 
to the effective action charge in the zero field limit, since 
the equation for Q is of Poisson type. Thus, e.g., the 
existence of dyons p(J will not cause a velocity shift in 
the weak field limit. 



APPENDIX A: 



Our aim is to evaluate the integral of eq.(Mq) 



1(h) 



dz 



-2hz 



; coth z — 1 1 



sinh 



coth ; 
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For this, we have to decompose it into simple parts which 
one can handle by standard methods of dimensional reg- 
ularization. Note that the integral is convergent, since 
the prescription h — > h — ie is implicitly understood. 

We begin with an integration by parts of the first term 
in square brackets with respect to the sinh 2 in the de- 
nominator. This leads to 



1(h) = I d:v' 
o 



-2hz 



2 1 

coth z — h coth z H - coth z 

2z z 



11 1 , 

t. coth z 

2z sinh 2 z 3 
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The last three terms of the expression in square brackets 
are already in a convenient shape. In the following, we 
thus consider only the remaining first two terms. The 
strategy is similar: we extract a term proportional to 
1/ sinh 2 z and integrate by parts. 



h(h) 



h J d:i-' 2hv 
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coth z 
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sinh 2 , 



h J dzc- 2hz [(2h+±)cothz+ (h+±)smh2z 



"(l+2^ + 1 ) COSh ^-(l+2^ + 1 )] 

(A3) 

Inserting Ji into eq. ( [A2] ) , we obtain the wanted types of 
integrals. Each of these can be integrated separately by 
introducing an extra factor of z e and rotating the contour 
onto the positive real axis. At the end, the 1/e-poles 
cancel and we arrive at the result given in eq.([53|) in the 
limit e — > 0. 



APPENDIX B: 

In this appendix, we want to expand the infinite sum 
in eq.(^) for small values of A := ^ which corresponds 
to a high temperature limit 



S(X) :=£(-!)" (An) 2 X 2 (An) 



(Bl) 



n=l 



Since the appearance of Bessel functions reflects the 
R 3 x S 1 topology which is the finite-temperature field 
theory space, the techniques described in the following 
are certainly useful for further finite-temperature appli- 
cations. 

The first step is to choose a representation of the mod- 
ified Bessel function that shows a simple dependence on 
the summation index |23|| 



^2 (An) 



-An cosh t 



cosh 2tdt, 
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Inserting eq.(|B2|) into eq.(Bl), leads us to 

°? oo 

S(X) = A 2 dt cosh2t^n 2 e- (i7r+Acosht) " 



(B3) 



By differentiating the geometric series X^^Lo 1™ = 
twice with respect to q, we find the result for the sum in 
eq-© 
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Inserting eq.(B4) into eq.(B3) and decomposing the 



cosh2i into 2 cosh t — 1, we get 
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where we have substituted p :— A coshi. With some care, 
the parameter integrals J\ and J 2 can now be expanded. 
We have to pay special attention to the process of taking 
the limit A — > for the J's and their derivatives. We can 
circumvent possible convergence problems at the lower 
bound by a repeated integration by parts of the square 
root terms. Using the short form 



1 (i- e -p) 
P (1 + e-f) 3 
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the non-vanishing coefficients of the expansion up to or- 
der C(A 5 ) can be expressed as 
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(BIO) 
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Finally, the Taylor expansion of eq.(Bl) reads 
5(A) = Ji(0) + (|Jf(0) + J 2 (0)) A 2 

+ (i Jf (0) + |4'(0)) A 4 + 0(A 6 ) 

= 1-^A 4 + 0(A 6 ), 

where the constant k\ is defined by 
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The constant &2 that appears in the calculation of the 
effective action charge for high temperature and strong 
fields is obtained by similar techniques. Its integral rep- 
resentation is (accidentally) equal to J%(0) 

J V (1 + e-f 3 
o 

= 0.213139199 408.... (B14) 



[23] I.S. Gradshteyn and I.M. Ryzhik, Tables of Integrals, Se- 
ries and Products, Academic Press (1965). 

[24] D. Bailin and A. Love, Introduction to Gauge Field The- 
ory, IOP Publishing Limited (1993). 

[25] G.K. Savvidy, Phys. Lett. B 71, 133 (1977); S.G. 
Matinyan and G.K. Savvidy, Nucl. Phys. B 134, 539 
(1978). 

[26] S.G. Kovalevich, P. Osland, Ya.M. Shnir and E.A. Tolka- 
chev, Phys. Rev. D 55, 5807 (1997). 



[1] S.L. Adler, Ann. Phys. (N.Y.) 67, 599 (1971). 
[2] E. Brezin and C. Itzykson, Phys. Rev. D 3, 618 (1971). 
[3] I.T. Drummond and S.J. Hathrell, Phys. Rev. D 22, 343 
(1980). 

[4] K. Scharnhorst, Phys. Lett. B 236, 354 (1990). 
[5] G. Barton, Phys. Lett. B 237, 559 1990. 
[6] R.D. Daniels and G.M. Shore, Nucl. Phys. B 425, 634 
(1994). 

[7] R.D. Daniels and G.M. Shore, Phys. Lett. B 367, 75 
(1996). 

[8] J.L. Latorre, P. Pascual and R. Tarrach, Nucl. Phys. B 

437, 60 (1995). 
[9] Wu-yang Tsai and T. Erber, Phys. Rev. D 12, 1132 

(1975). 

[10] G.M. Shore, Nucl. Phys. B 460, 379 (1996). 

[11] G. Barton and K. Scharnhorst, J. Phys A 26, 2037 

(1993) . 

[12] W. Dittrich and M. Reuter, Effective Lagrangians in 
Quantum Electrodynamics, Lecture Notes in Physics 220, 
Springer- Verlag (1985). 

[13] J. Schwinger, Phys. Rev. 82, 664 (1951). 

[14] D. Bakalov et at, Nucl. Phys. B 35 (Proc. Suppl), 180 

(1994) . 

[15] D. Bakalov et al, INFN/AE-94/27 unpublished (1994). 
[16] E. Iacopini and E. Zavattini, Phys. Lett. B 85, 151 
(1979). 

[17] V.I. Ritus, JETP 42, 774 (1975); Issues in Intense-Field 
Quantum Electrodynamics, Proc. Lebedev Phys. Inst. 
168, New York (1987); 

for a modern treatment employing the string-inspired 
world line formalism, cf. M. Reuter, M.G. Schmidt and 
C. Schubert, Ann. Phys. (NY) 259 313 (1997); B. Kors 
and M.G. Schmidt, [hep-th/9803144| unpublished (1998). 

[18] W. Dittrich and H. Gies, unpublished; W. Dittrich, QED 
in Constant External Magnetic Fields, Lecture Notes 
(WS 1997/98), Univ. Tubingen (1997). 

[19] L. Bendersky, Acta Mathematica (Uppsala) 61, 263 
(1933). 

[20] L.S. Brown and G.J. Maclay, Phys. Rev. 184, 1272 
(1969). 

[21] X. Kong and F. R avndal, Phys. Rev. Lett. 79,545 (1997); 

[hep-ph/9803216| unpublished (1998). 
[22] W. Dittrich, Phys. Rev. D 19, 2385 (1979). 



13 



